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Abstract 

We propose a parametrization of the lepton mixing matrix in terms of an 
expansion in powers of the deviations of the reactor, solar and atmospheric mix- 
ing angles from their tri-bimaximal values. We show that unitarity triangles and 
neutrino oscillation formulae have a very compact form when expressed in this 
parametrization, resulting in considerable simplifications when dealing with neu- 
trino phenomenology. The parametrization, which is completely general, should 
help to establish possible relations between the deviations of the reactor, solar 
and atmospheric mixing angles from their tri-bimaximal values, and hence enable 
models which predict such relations to be more directly compared to experiment. 



E-mail: sfk@hep.phys.soton.ac.uk 



Over the last decade neutrino physics has undergone a revolution with the mea- 
surement of neutrino mass and lepton mixing from a variety of solar, atmospheric and 
terrestrial neutrino oscillation experiments [l]. Lepton mixing is described by the 3x3 
matrix [2] 
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The Particle Data Group (PDG) parameterization of the lepton mixing matrix (see e.g. 
[S])is: 

(C12C13 S12C13 sise""^ \ 

-C23S12 - Si3S23Ci2e"^ C23C12 - Si3S23Sl2e"' S23C13 -P , (2) 
S23S12 - Si3C23Ci2e'^ -S23C12 - Si3C23Si2e"' C23C13 / 

where S13 = sin 6^13, C13 = cos6'i3 with 613 being the reactor angle, S12 = sin6'i2, 
C12 = cos ^12 with 6*12 being the solar angle, S23 = sin 6^23, C23 = cos 6*23 with 6*23 be- 
ing the atmospheric angle, 6 is the (Dirac) CP violating phase which is in principle 

measurable in neutrino oscillation experiments, and P = diag(e' 2 2 ,0) contains ad- 
ditional (Majorana) CP violating phases ai,a2- Current data is consistent with the 
tri-bimaximal mixing (TBM) formal [1] 
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Many models can account for TBM lepton mixing [SI El [TJ El El [lOl [H]- However there 
is no convincing reason for TBM to be exact, and in the future deviations from it are 
expected to be observed. With this in mind it is clearly useful to develop a parametriza- 
tion of the lepton mixing matrix in which such deviations are manifest, and in which 
the predictions of models for deviations from tri-bimaximal mixing can naturally be 
expressed. Such a parametrization must be model independent, and completely general 
so that it can be used by experimentalists and phenomenologists in performing analyses 
of neutrino experiments. It must also be sufficiently simple to be useful and yet accurate 
enough to be reliable. 

In this paper we discuss a parametrization of the lepton mixing matrix which pos- 
sesses all of the above desirable features. The parametrization exploits the empirical 
observed closeness of lepton mixing to the TBM form, and is analagous to the Wolfen- 
stein parametrization of quark mixing [T2] . Just as the Wolfenstein parametrization is an 



^Sometimes an alternative phase convention is chosen in which the third row of ?7mns has its signs 
reversed. 
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expansion about the unit matrix, so the present parametrization is an expansion about 
the tri-bimaximal matrix. Unhke the Wolfenstein parametrization, we introduce three 
small parameters parametrizing the deviations of the reactor, solar and atmospheric 
angles from their tri-bimaximal values. The expansion works since all three parameters 
are empirically small, having magnitude of order the Wolfenstein parameter A ~ 0.227 
or less. A related proposal to expand the lepton mixing matrix elements about the 
tri-bimaximal matrix elements, using a different parametrization from that introduced 
here, was discussed in p3]. Other related proposals to parametrize the lepton mixing 
matrix have been considered in [TH [151 HBl [13, [IE] El 
Without loss of generality we define 

r 1 1 

si3 = si2 = -^(l + s), S23 = -j={l + a), (4) 

where we have introduced the three real parameters r, s, a to describe the deviations of 
the reactor, solar and atmospheric angles from their tri-bimaximal values. Global fits 
of the conventional mixing angles [19] can be translated into the 2a ranges H 

< r < 0.22, -0.11 <s< 0.04, -0.12 <a< 0.13. (5) 

The empirical smallness of these parameters suggests that we consider an expansion of 
the lepton mixing matrix in powers of r, s, a about the tri-bimaximal form. To first 
orderly in r, s, a the lepton mixing matrix can be written 



U 



-^(1 + s - a + re'^) -1.(1 - - a - \Te'') ^(1 + a) 

V 75(1 + ^ + « - ^^'') -73(1 - + « + 75(1 - / 



P. (6) 



As in the Wolfenstein parametrization, the above parametrization of the lepton mixing 
matrix avoids the introduction of mixing angles, instead dealing directly with elements 
of the mixing matrix. Accordingly the parametrization results in considerable simplifi- 
cations when dealing with neutrino phenomenology. For example, the complex elements 
of the quark mixing matrix can be visualized using unitarity triangles [20], which, when 
normalized, only depend on two parameters. The same proves to be true when using the 
above parametrization of the lepton mixing matrix. The sides of the unitarity triangles 
enter into the neutrino oscillation formulae, and consequently these are also considerably 



am grateful to S. Parke, Z.Z. Zing and P. Harrison for informing me about their work [T4 l [TSl fl7] . 
^Note that r must be positive definite, while s, a can take either sign. Indeed there is a preference 
for s to be negative. 

"^The second order corrections are expected to be very small, of order one per cent or less, depending 
on the (presently constrained but undetermined) values of r, s, a. Throughout the main text where 
results are presented to first order in r, s, a, the second order corrections are given in Appendix O 
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Figure 1: The V2.v^ unitarity triangle. The angle 7 is equal to the CP phase 5 to first order. The 
unknown Majorana phases just rotate the triangle in the complex plane. The rescaled triangle is 
oriented as shown with the opening angles unchanged, the horizontal side having unit length, and the 
shortest side having length r to first order. Currently < r < 0.22 at 2cr, and the opening angles a, (3 
and 7 are all undetermined. 



simplified by the new parametrization. In the remainder of the paper we shall discuss 
unitarity triangles and neutrino oscillation formulae using the above parametrization. 
CP violation is described by the Jarlskog [2T] invariant which to leading order is 

J^Uind. (7) 

D 

Leptonic unitarity triangles [22] may be constructed using the orthogonality of different 
pairs of columns or rows of the mixing matrix. Only the opening angles, side lengths and 
areas of the triangles have physical significance. For example the area of each unitarity 
triangle is A = ^\J\ and CP violation implies that the longest side of each unitarity 
triangle is smaller than the sum of the other two. Current solar, reactor and atmospheric 
experiments directly constrain the elements Ue2, f^es and which have a particularly 
simple parametrization in EqEl The most important unitarity triangles should therefore 
include all of the elements Ue2, t^es and ?7^3. There are two such unitarity triangles, the 
Z/2.Z/3 one [in] corresponding to the orthogonality of the second and third column, and 
the i^e-^fi one [23] corresponding to the orthogonality of the first and second row. Each 
of them has a simple expression in terms of the new parametrization, as we now discuss. 
The z/2.^'3 triangle in FigH] corresponds to the unitarity relation 

Ue2U:, + U,2U;, + Ur2U:, = 0. (8) 

To first order the sides of this unitarity triangle are given by 
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Ss = Ur2U:,^-^{l-'- + '-e^'). (9) 
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Figure 2: The Ve-i^fi unitarity triangle. The angle 7' is equal to the CP phase 6 to first order. The 
unknown Majorana phases cancel. The rescaled triangle is oriented as shown with the opening angles 
unchanged, the horizontal side having unit length, and the shortest side having length |r to first order. 
Currently < r < 0.22 at 2(t and the opening angles a' , (3' and 7' are all undetermined. 

Clearly + ^2 + ^3 = to first order. The invariant J is 

J = Im{SiS;) = ImiSsSl) = Im{S2S;) (10) 

which yields EqJTl To first order the sides of this triangle are only sensitive to the solar 
and reactor parameters s and r and the phase 6, with the atmospheric parameter a only 
appearing at second order. One may rescale the sides by 5*3 
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To first order the rescaled triangle is only sensitive to the reactor parameter r and the 
phase S, which is the anticipated result. To second order the solar parameter s (but not 
the atmospheric parameter a) appears. 

The other unitarity triangle of interest is I'e-'^ii in Fig J21 corresponding to the unitarity 
relation 



(12) 



To first order the sides of this unitarity triangle are given by 

T3 = u,su:. 
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Clearly Ti + T2 + T3 = to first order. The invariant J is 



J = Iminr;) = Im{TiT;) = Im{T2T*) (14) 

which again yields EqjTl Unlike the previous case, the sides of this triangle are sensitive 
to the atmospheric parameter a at first order. One may rescale the sides by Ti 

T[ = 1 

^/ _ Ufi2U*2 _ _i , 3 i6 

As in the previous case, to first order the rescaled triangle is only sensitive to the reactor 
parameter r and the phase S, which is the anticipated result. To second order the solar 
parameter s and the atmospheric parameter a appear. 

We now turn to the application of the parametrization in EqHJto neutrino oscilla- 
tions. Let us denote by Pajj = P^i^a ^p) the probability of transition from a neutrino 
fiavour a to a neutrino fiavour (3. Then expanding to second order in the parame- 
ters r, s, a and A21, where it is assumed that A21 ^ 1 as in [24j, we find considerably 
simplified vacuum oscillation [f| probabilities. 

The electron anti-neutrino disappearance probability relevant for a reactor experi- 
ment |25] is given to second order in r, s, a and A21 as 

Pee = 1 - sin^ A31 - ^A^, (16) 

where Ajj = 1.27Am^jL/E with L the oscillation length in km, E the beam energy 
in GeV, and Amfj = mf — in eV^. Note that this disappearance probability is 
independent of the solar and atmospheric parameters s, a, as well as the phase 6, to this 
order. 

The electron neutrino appearance probability relevant for a forthcoming long baseline 
muon neutrino beam experiment [26j is given to second order in r, s, a and A21 as 

4 4 

P^e = r^sin^Asi + -A^^ + -rA2isinA3icos(A3i + 5). (17) 

It is also independent of the solar and atmospheric parameters s, a and only depends on 
the reactor parameter r and the phase 6 to this order. The reason is that each of the 
terms is second order in the parameters r, A21, so any deviations from tri-bimaximal solar 



^Similar considerations apply to oscillations in the presence of matter as discussed in Appendix [Bl 
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or atmospheric mixing only appear at third order. The muon neutrino disappearance 
probabihty is given to second order in r, s, a and A21 as 

P^^ = l-(l-4a2)sin2A3i-^(l + 3cos2A3i)A2, 
2 

+ -(1 -s-rcos5)A2isin2A3i. (18) 

Muon neutrino disappearance is clearly sensitive to deviations from tri-bimaximal mix- 
ing, since all three parameters r, s, a and the phase 6 appear. For example the prospects 
for measuring deviations from maximal atmospheric mixing in the next generation of 
long baseline muon neutrino beam experiments has recently been discussed [27]. Simi- 
larly the tau neutrino appearance probability is given to second order in r, s, a and A21 
as 

P^r = (l-4a2-r2)sin2A3i-^(l-3cos2A3i)A^i 

- ^(l-s)A2isin2A3i + -rA2isin2A3isin(5. (19) 

We emphasize that the parametrization discussed here is completely general and 
is not based on the ansatz of tri-bimaximal mixing, any more than the Wolfenstein 
parametrization [12j is based on the ansatz that the quark mixing matrix is equal to the 
unit matrix. Just as the Wolfenstein parametrization is an expansion about the unit 
matrix, so this parametrization is an expansion about the tri-bimaximal matrix. Unlike 
the Wolfenstein parametrization, there are three small parameters r, s, a parametrizing 
the reactor, solar and atmospheric deviations from tri-bimaximal mixing. The expansion 
works since the deviations from tri-bimaximal mixing are empirically small parameters 
with r, s, a all having magnitude of order the Wolfenstein parameter A ~ 0.227 or less. 
Indeed these parameters are sufficiently small that the first order approximation is accu- 
rate enough for many purposes, resulting in quite a simple looking lepton mixing matrix 
in EqEl for example. Unitarity triangles and neutrino oscillation formulae also have a 
very simple form when expressed in this parametrization. 

The three parameters r, s, a are not determined at the present time, and it is even 
possible that one or more of them (possibly all of them) are zero, although this seems 
a priori unlikely. However, as mentioned, many speculations appear in the literature 
as to the origin and nature of tri-bimaximal mixing and the deviations from it, and 
these speculations naturally find expression in this parametrization. For example certain 
classes of unified fiavour models |5j predict a sum rule which relates s to r and 6, namely 
s ~ rcosS, where r A/3 and a = O(A^). Alternatively it has been suggested |16j that 
trimaximal solar mixing is exact, s = 0, with a ~ —^rcos6 and r unspecified. Clearly 
an important goal of the next generation of neutrino experiments must be to show that 
the parameters r, s, a differ from zero. Subsequent high precision neutrino experiments 
will then be required to accurately measure the values of the parameters 5, Ci, clS well 



6 



as (5, to investigate their possible relationships to each other and to the Wolfenstein 
parameter A. 



Acknowledgments 

We would like to thank J. Flynn for useful discussions. We acknowledge partial support 
from the following grants: PPARC Rolling Grant PPA/G/S/2003/00096; EU Network 
MRTN-CT-2004-503369; EU ILIAS RII3-CT-2004-506222. 

Appendix 

A Second order corrections 

In this appendix we list the second order corrections to all the results given in the main 
text. The second order corrections to the first order MNS matrix elements in Eqjnjare 




v/3' 2 2 2 8 




V3'2 2 2 8 




(20) 



The second order correction to the Jarlskog CP invariant in EqJT] is 



A J ^ — sin 6. 



rs 



(21) 
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The second order corrections to the unsealed sides of the 1^2 •^'3 unitarity triangle in 
EqJH] are 

AS2 ^ -^C- + 2a^ + ls^ + are'' + lsre'') 
V6 4 8 2 

AS, ^ l=(^!l + 2a' + ^s' + are''-^sre''). (22) 

The second order corrections to the normalized sides of the i'2-^3 unitarity triangle in 
EqJTTl are 

AS[ ^ !le''^-^sre'^ 
^ 2 2 

AS', ^ -i^^^ + Lre^^ 

AS's = 0. (23) 

The second order corrections to the unsealed sides of the i^e-^fj, unitarity triangle in 
EqJTSl are 

« ^ 7 r, a"^ sa sr -c or 

ATi ^ — + —s^ + — + — + —e'^ e'^ 

^ 12 24 3 6 3 3 

AT2 ^ e'^ e*^ 

12 24 3 6 3 6 

ATg ^ ye^^. (24) 

The second order corrections to the normalized sides of the z/g.z/^ unitarity triangle in 
EqJTSl are 

at; = 

AT^ ^ -^re'\2re'^ + s - Aa) 

ATg ^ ^re'\2re'^ + s-Aa). (25) 

B Neutrino oscillations in matter 

In this appendix we present the complete formulae for neutrino oscillations in the pres- 
ence of matter of constant density to second order in the quantities r, s, a and A21, where 



8 



it is assumed that A21 ^ 1 as in t24j. Following [24j we write A = A31, a = =^ and 
A = ^ where V is the potential expressed in units of eV as 

V ^ 7.56 X Ye (26) 

where p is the matter density of the Earth in units of g/cm^ and ^ 0.5 is the number 
of electrons per nucleon in the Earth. The constant density approximation is good 
when the neutrino beam only passes through the Earth's crust where p ~ 3 g/cm^ or 
the Earth's mantle where p ~ 4.5 g/cm^. 

The complete set of neutrino oscillation probabilities for electron neutrino or muon 
neutrino beams in the presence of matter of constant density to second order in the 
parameters r, s, a and a are 

8 gSinMA 2sin^(^-l)A 
4 ,sinMA ,sin2(A-l)A 

9 ^ {A-iy 

4 ,^ sinAAsin(v4- 1)A 
+ -racos(A-5) ^- (28) 

4 ,sinMA ,sin^(A-l)A 
Per = -a n^ + r^ 
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9 A^ (A-i; 

4 ^,sinAAsin(yl- 1)A 
- 3™co.(A-*)^ i^-^ (29) 

4 ,sinMA ,sin^(v4-l)A 

9 (A -1)2 

4 siny4Asin(A- 1)A , , 

+ -mcos(A + 5)^^ p^T)^- 
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1 - (1 - Aa^) sin^ A + - (1 - s)aAsin2A 
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4 ^sinMA 4 

2 --a2A2cos2A 
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4 2I , . A 

-a — smA 

9 A 
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sin A A , , X . A 

cos(A - 1)A sin2A 
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— r 



,sm^{A - 1)A 



smAcosAA — ^ ^ Asin2A 

A - 1 V - 1) 2 

4 ^ ^ sin AA sin(A - 1)A 
-ra coso cos A — — . 

3 A {A-\) 



(31) 



2 4 

P^^ = (1 -4a^)sin2 A - -(1 - s)aAsin2A + -a2A2cos2A 

3 9 

4,1/ sin A A , ^ , , A 

- -a^— sinA — cos A - 1)A sin 2 A 

9 A V A ^ ' 2 

1 2 . . . sinfA - 1) A A ^ 

I smAcosAA — ^. ^ —Asm 2 A 



+ 



A-1 



(A-l) 
sin AA sin(A - 1)A 



+ -msin^sinA ^ _ 



(32) 
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